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arise through tree-level exchanges of heavy particles, or through strong-interaction 
loop diagrams. The Wilson coefficients of the corresponding electroweak penguin 
operators are calculated in a large class of New Physics models, and in many cases 
are found not to be suppressed with respect to the QCD penguin coefficients. Sev- 
eral tests for these effects using observables in i? ± — >• nK decays are discussed, and 
nontrivial bounds on the couplings of the various New Physics models are derived. 
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1. Introduction 

The study of rare decay processes is an important tool in testing the fundamental 
interactions among elementary particles, exploring the origin of CP violation, and 
searching for New Physics beyond the Standard Model (SM). Such processes have 
been explored in great detail, both theoretically and in experimental searches, in the 
weak interactions of kaons and B mesons, as well as in K-K, D-D, and B-B mixing. 
Two prominent examples are the evidence for mixing-induced CP violating in the 
decay B — > J/ipK s reported by the CDF Collaboration [1], and the observation of 
direct CP violation in the decays K — > tttc [2, 3], which has recently been confirmed 
by the KTeV and NA48 Collaborations [4, 5]. 

Flavor-changing neutral current (FCNC) processes, which are forbidden at tree 
level in the SM, are especially sensitive to any new source of flavor-violating in- 
teractions. Already, the absence of experimental signals for new FCNC couplings 
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puts stringent bounds on the parameters of many extensions of the SM such as 
Supersymmetry (SUSY) [6]. So far, FCNC processes have been explored mainly 
in particle-antiparticle mixing and in "semihadronic" weak decays, which permit a 
clean theoretical description. In the kaon system, examples of the latter type are the 
decays K — > 7r and K — > 7r vv. In the B system, the decays that have received 
the most attention are B — > X s 7, B — > X s / + Z~ and 5 — > X s i/P, where X s can be 
any final state, exclusive or inclusive, containing a strange quark [7]. 

In the present paper, we explore in detail how New Physics could affect purely 
hadronic decays such as B — > ttK, which are sensitive to isospin- or, more generally, 
SU(3) flavor-violating interactions. In the SM, the main contributions to the de- 
cay amplitudes for these processes come from the penguin-induced FCNC transition 
b — > sqq, which by far exceeds a small, Cabibbo-suppressed b — > uus contribution 
from lU-boson exchange. Because of a fortunate interplay of isospin, Fierz and flavor 
symmetries, the theoretical description of the charged decays B^ — > nK in the SM 
is clean despite the fact that these are exclusive nonleptonic decays [8, 9]. Isospin 
violation arises through the small charged- current contribution and through the elec- 
troweak penguin operators in the low-energy effective weak Hamiltonian [10]. In the 
SM, these operators are induced by penguin and box diagrams involving the ex- 
change of weakly interacting W and Z bosons, or of a photon. Here we point out 
that in a large class of New Physics models such effects can be mediated by "trojan" 
electroweak penguins, which are neither (pure) penguins nor of electroweak origin. 
Nevertheless, at low energies their effects are parameterized by an extension of the 
usual basis of electroweak penguin operators. We will explore examples where trojan 
penguins are induced by tree- level couplings, e.g., models with an extra Z' boson, or 
pure strong- interaction processes, e.g., gluino box diagrams in SUSY models. 

In the present work we calculate the Wilson coefficients of the hadronic elec- 
troweak penguin operators in the effective weak Hamiltonian in an extended operator 
basis and for a large class of New Physics models. We then explore the phenomeno- 
logical consequences of these new, isospin- violating contributions for weak-interaction 
observables. This serves two purposes: first, it allows us to derive bounds on New 
Physics parameters, which in some cases improve upon existing bounds derived from 
other processes; secondly, it shows which observables may be interesting to look at 
as far as searches for New Physics are concerned. We shall address both issues in 
detail for the particular case of B ± — > nK decays, which in the context of the SM are 
of prime importance in determining the weak phase 7 = arg(U M * 6 ) of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [8, 9, 11, 12, 13, 14]. We stress, however, that 
trojan penguins could be important in a much wider class of processes. In particular, 
they may be responsible for a large contribution to the quantity e' / e measuring direct 
CP violation in K — > decays [15, 16]. The importance of B^ — > ttK decays in the 
search for New Physics has been emphasized in [9, 17], and some specific scenarios 
containing new isospin- violating contributions have been explored in [18, 19]. 
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In Section 2, we discuss the effective weak Hamiltonian relevant to hadronic B 
decays. In the presence of generic New Physics contributions, the basis of penguin 
operators has to be extended from the standard one in several aspects. We discuss the 
structure of the new operators and their scaling properties under a renormalization- 
group transformation from a high scale down to low energies. The theory of the 
rare hadronic decays — > ixK is discussed in Section 3, where we indicate how 
various observables in these decays can be used to test for physics beyond the SM. 
We derive mo del- independent bounds on a ratio of CP-averaged B ± — > nK 
branching ratios in the presence of New Physics, and show that the value 7^ of the 
weak phase extracted in B — > %K decays is extremely sensitive to isospin-violating 
New Physics contributions. Because the theoretical analysis has only small hadronic 
uncertainties, potential New Physics effects can be detected even if they are 50% 
smaller than the isospin-violating contributions present in the SM. In Sections 4 and 
5, we calculate the Wilson coefficients of the b — * sqq penguin operators in a large 
class of extensions of the SM, including models with tree-level FCNC couplings of the 
Z boson, extended gauge models, multi-Higgs models, and SUSY models with and 
without R-parity conservation. In each case, we explore which region of parameter 
space can be probed by measuring certain B ± — > nK observables, and how big a 
departure from the SM predictions one can expect under realistic circumstances. 
Section 6 contains a summary of our results and the conclusions. 



2. Effective Hamiltonian for hadronic FCNC processes 

The effective weak Hamiltonians relevant to rare hadronic B decays based on the 
quark transitions b — > sqq or b — > dqq have been discussed extensively in the litera- 
ture. Here we consider only the first case, adopting the notations of [10]. A similar 
discussion (with obvious replacements of indices) would apply to the other case. In 
the SM, the result can be written in the compact form 



n ^ = % E \ [ E Q q M + E CiMQiM + cMQM 

q=u,c L i=l,2 i=3.. .10 



, (2-1) 



where X q = V* b V qs are combinations of CKM matrix elements obeying the unitarity 
relation \ u + \ c + \ t = 0, and Qi are local operators containing quark and gluon 
fields. Specifically, 

Q\ = (hqp)v-A (qpS a ) v _ A , Q\ = (b a q a ) v _ A (g /3 s /3 )y_ A , (2.2) 

summed over color indices a and /3, are the usual current-current operators induced 
by VU-boson exchange, 

Q3 = (b a s a ) v _ A ^{qpq^v-A, 
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Qi — {b a S/3) V -A ^(QpQcdv-A, 

q 

Q 5 = (b a s a )v-A ^2(q/3Q/s)v+a , 
q 

Q% = (b a Sp) V -A ^2{qpq a )v+A, (2.3) 
q 

summed over the light flavors q = u,d, s, c, b, are referred to as QCD penguin oper- 
ators, and 

3 

Qi = - (b a s a ) v ^ A e q (qpvph+A , 

1 q 
3 

Qs = - (b a sp)v-A e q iw°)v+A , 
A q 

3 

Q9 = - (b a s a ) v „ A J2 e q {wp)v-a , 

1 q 
3 

Qio = -z (b a sp) V -A e q i.Wo)v~A , (2.4) 
1 q 

with e q denoting the electric charges of the quarks, are called electroweak penguin 
operators. The notation (qiq2)v±A implies <7i7 M (l ±75)^2- The terminology of QCD 
and electroweak penguins is slightly misleading insofar as the Wilson coefficients 
of the QCD penguin operators also receive small contributions from electroweak 
penguin and box diagrams. However, in the SM there are no strong-interaction 
contributions to the coefficients of the electroweak penguin operators. The opera- 
tor Q 8g = 2 ^-ba^ u (l — 75)G M ^s in (2.1) is the chromomagnetic dipole operator. 
The analogous electromagnetic dipole operator and semileptonic operators contain- 
ing products of a quark current with a lepton current can be safely discarded from 
the effective Hamiltonian for hadronic B decays. 

In general, physics beyond the SM can induce a much larger set of penguin 
operators, and it is therefore unavoidable for our purposes to generalize the standard 
nomenclature reviewed above. However, for all the models we explore below it is 
sufficient to consider products of vector and/or axial vector currents only. We define 
a basis of such operators by 

Of = {b a s a ) V -A (qpq^v+A , 0\ = (b a sp) V -A {.qpq a )v+A , 

01 = {b a s a )v^ A (qpqp) V -A , 0\ = {b a Sp) v _ A (qpq a )v-A , 

Oi = (b a q a )v-A (qpSp) v+A , 0% = {b a q p ) v ^ A (qpS a ) v+A , (2.5) 

and denote their Wilson coefficients by cf. We implicitly assume a regularization 
scheme that preserves Fierz identities. In a general model, we also need operators 
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of opposite chirality compared to the ones shown above, i.e., with V — A <-> V + A 
everywhere. We denote these operators by Of and their coefficients by cf . Thus, the 
most general penguin Hamiltonian considered in this paper takes the form 

^Peng = ^ E E [*) + • (2-6) 

To this, one has to add the current-current operators in (2.2) and the chromo- 
magnetic dipole operator. It is implicitly understood that for the cases where q = s 
or b the operators Of and Of as well as Of and Of are omitted from the list of 
operators, because they are Fierz-equivalent to the remaining ones. 

The QCD and electroweak penguin operators present in the effective weak Hamil- 
tonian of the SM in(2.1) are linear combinations of the four operators 0\ A . However, 
in a general model there may be additional penguin operators built out of Of, Of, and 
of the opposite-chirality operators Of. Also, there may be operators which cannot 
be represented as linear combinations of QCD and electroweak penguin operators as 
shown in (2.3) and (2.4), because in a general model the flavor symmetry among up- 
or down-type quarks can be violated. An example would be operators with flavor 
content sb(dd — ss), which are absent in the SM. Because such operators are of no 
relevance to our discussion they will not be explored any further here. For complete- 
ness, we show how the Wilson coefficients O 3 ...i of the SM penguin operators in (2.1) 
can be expressed in terms of the coefficients cf . Defining the linear combinations 

C Q^1±M, cF w = cr-cf, (2.7) 



we have 



— ^<^3,4 — C3 4 , — A t 5i6 — C 12 

-lx t a = cf w 6 ; i = 7,. ..,10. (2.8) 

The isospin- violating effects induced by the coefficients cf w and cf w are the main 
focus of this paper. In the SM, the matching conditions for the corresponding elec- 
troweak penguin coefficients at the weak scale, and to leading order in perturbation 
theory, are C 8 (mw) = Cio( m w) — 0, Cj(m w ) m 0, and [20] 

C s (m w ) » ^ (l + , (2.9) 

12-7T sm B w \ x t — 1 J 

where Xt = {m t /mw) 2 - For simplicity, we show only the large electroweak contri- 
bution to Cg(mw) and omit a common, renormalization-scheme dependent electro- 
magnetic contribution to C 7 (m w ) and Cg(mw), which is negligible compared with 
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the contribution in (2.9). In this approximation, the SM matching conditions for the 
electroweak penguin coefficients in our basis read 



c^ SM (m w ) = f ( 1 + ^ ) - -5.5 x 10- , 

87T sin ^ y x t — 1 J 

EW,SM/ \ ~EW,SM/ \ n /a-,M 

c^ 3 (mif) = Ci {m w ) = 0. (2.10) 

For the numerical estimate we have used m t = m t {m t ) = 170 GeV, a = 1/129, and 
A t = V£V U = -0.04. 

In Sections 4 and 5, we will calculate the Wilson coefficients cf at a high scale, 
which for simplicity will be identified with the electroweak scale. In phenomenological 
applications, however, one usually prefers working with coefficients renormalized at a 
low scale of order m b . The two sets of coefficients are connected by a renormalization- 
group transformation [10]. Here we discuss the QCD evolution of the electroweak 
penguin operators at the leading logarithmic order, neglecting next-to-leading cor- 
rections as well as QED corrections to the anomalous dimensions of the operators. 
Under QCD evolution, flavor-nonsinglet combinations of penguin operators mix into 
flavor- singlet combinations through diagrams in which two light quarks annihilate 
into a gluon, which then fragments into a pair of light quarks. However, there is no 
mixing of flavor- singlet operators into flavor-nonsinglet ones. If we restrict ourselves 
to flavor-nonsinglet combinations of the operators Of, each pair in the three lines 
in (2.5) obeys a separate matrix evolution equation. It follows that the coefficients 
cf w mix pairwise under renormalization. For each pair of coefficients associated with 
(V =f A) ® (V ± A) operators we obtain 



k-24/23 _ „-3/23 

cf W M = ^ c™(m w ) + ^ c™(m w ) , (2.11) 



and similarly for (cf w ,cf w ), where k = a s (ii)/a s (mw)- For a pair of coefficients 
associated with (V =F A) <g> (V =F A) operators we obtain instead 

k-12/23 , K -6/23 K 12/23 _ ,.-6/23 

cf W M = ^ <T{m w ) - ^ c™(m w ) , 

„ 12/23 _ K -6/23 k-12/23 , K -6/23 

c 4 EW (//) = ^ c™(m w ) + K - ±H c 4 EW (™vk) . (2.12) 

The coefficients cf scale in the same way as the cf. Since our main focus is on 
electroweak penguins and their generalizations beyond the SM, we will not discuss 
the more complicated evolution equations for the coefficients cf themselves, which 
can however readily be deduced from the literature [10]. 
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3. Searching for New Physics with B 



nK decays 



In close correspondence with the different types of operators in the effective weak 
Hamiltonian, one distinguishes three classes of flavor topologies relevant to B — > ttK 
decays, referred to as trees, QCD penguins and electroweak penguins. In the SM, 
the weak couplings associated with these topologies are known. From the measured 
branching ratios for the various B — > nK decay modes it follows that the QCD 
penguins dominate the decay amplitudes [21], whereas trees and electroweak pen- 
guins are subleading and of a similar strength [22]. The theoretical description of 
the two charged modes B ± — > n ± K° and B* 1 — > n^K^ exploits the fact that the 
amplitudes for these processes differ in a pure isospin amplitude A 3 / 2 , defined as 
the matrix element of the isovector part of the effective Hamiltonian between a B 
meson and the nK isospin eigenstate with I — |. In the SM the parameters of this 
amplitude are determined, up to an overall strong-interaction phase 0, in the limit 
of SU(3) flavor symmetry [8]. SU(3)-breaking corrections can be calculated in the 
factorization approximation [23], so that theoretical uncertainties enter only at the 
level of nonfactorizable SU(3)-breaking corrections to a subleading decay amplitude. 
Moreover, it has recently been shown that even these nonfactorizable corrections can 
be calculated in a model-independent way up to terms that are power suppressed in 
A/rrib and vanish in the heavy-quark limit [24]. 



3.1 General parametrization of the decay amplitudes 

In the presence of New Physics, the analysis of B ± — > nK decays becomes more 
complicated. A convenient and completely general parametrization of the two decay 
amplitudes is 

A{B + -> n + K°) = P (1 - ipe**') , 
-V2A(B + -> n°K+) =P[l- ipe i(pp - e 3/2 e^(e* 7 -oe* - ibe i4 ">)] , (3.1) 

where P is the dominant penguin amplitude defined as the sum of all CP-conserving 
terms in the B^ — > n ± K° decay amplitudes, £3/2, p, a, b are real hadronic parameters, 
and 0, (f> p , (p a , (f>b are strong-interaction phases. The weak phase 7 = arg(V^* fe ) 
and the terms ip and ib change sign under a CP transformation, whereas all other 
parameters stay invariant. The terms proportional to e 3 / 2 in (3.1) parameterize the 
isospin amplitude A 3 / 2 . The contribution proportional to e n comes from the matrix 
elements of the current-current operators Qi and Q 2 in the effective Hamiltonian, 
which mediate the tree process b — > uus. The quantities a and b parameterize the 
effects of electroweak penguins. It is crucial that only isospin-violating terms can 
contribute to the amplitude A 3 / 2 . All isospin-conserving contributions reside in P 
and p. 
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Let us discuss the various terms entering the decay amplitudes in detail. The 
parameter e 3 / 2 characterizes the relative strength of tree and QCD penguin contri- 
butions. Information about it can be derived by using SU(3) flavor symmetry to 
relate the tree contribution to the isospin amplitude A3/2 to the corresponding con- 
tribution in the decay B + — > 7r + 7r°. Since the final state 7r + 7r° has isospin 1 = 2 
(because of Bose symmetry), the amplitude for the latter process does not receive 
any contribution from QCD penguins. Moreover, in the SM electroweak penguins 
in b — > dqq transitions are negligible, and thus only the tree topology contributes to 
the B + — > 7r + 7r decay amplitude. In our analysis we make the plausible assumption 
that potential New Physics contributions to this amplitude can be neglected. 1 Even 
if new electroweak penguin effects would be of comparable strength in b — > sqq and 
b — > dqq transitions, the latter would have to compete with a Cabibbo-enhanced tree 
amplitude in order to be significant in B + — > 7r + 7r° decays. We then find that [8, 9] 



£3/2 



£3/2 



V2R; 



SU(3) 



V u . 



V, 



ud 



B(B+ 



0^ 



B(B- 



7T-7T ) 



B(B+ -> n+K°) + B(B~ -> ir~K°) 



1/2 



(3-2) 



SU(3)-breaking corrections are described by the factor Rsv(3) — 1-22 ± 0.05, which 
can be calculated in a model-independent way using the QCD factorization theorem 
of [24]. The quoted error is an estimate of the theoretical uncertainty due to un- 
controllable corrections of O(j^^f-). Using preliminary data reported by the CLEO 
Collaboration [25] to evaluate the ratio of branching ratios in (3.2), we obtain 



£3/2 = 0.21 ± 0.06 cxp ± 0.01 th 



(3.3) 



With a better measurement of the branching ratios the uncertainty in £3/2 will be 
reduced significantly. 

The parameter p in (3.1) parameterizes the sum of all CP- violating contributions 
to the B + — > ir + K° decay amplitude. In the presence of New Physics, those contri- 
butions could arise from QCD as well as electroweak penguin operators. Note that 
the CP-conserving part of such terms is absorbed, by definition, into the quantity 
P. We will not attempt a theoretical calculation of this quantity (which is difficult 
even in the SM) and only consider observables that are independent of P. In the 
SM, p ~ e a sin7 [9] describes a small contribution induced by final-state rescattering 
from tree or annihilation diagrams [26, 27, 28, 29, 30, 31]. In the heavy-quark limit, 
the parameter e a can be calculated and is found to be of order —2% [24]. 

Finally, in the SM the parameter b vanishes, while 



ae* = Sew = (0.64 ± 0.09) x 



0.085 



\V ub /V c 



cb 



(3.4) 



1 If this were not the case, the New Physics impact on e 3 / 2 would provide us with another handle 
on non-standard isospin- violating effects. 
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is calculable in terms of fundamental parameters [8, 9, 32]. Up to some small SU(3)- 
breaking corrections, <5 E \y is given by the Wilson coefficient c^ w (m w ) in (2.10) di- 
vided by — |A M |. There are no additional hadronic uncertainties in this estimate in 
the SM. In particular, the strong-interaction phase <p a is bounded to be less than a 
few degrees and can be neglected for all practical purposes [9]. In a general model, 
the parameters a and b depend on the values of the penguin coefficients cf w and 
cf w as well as on the hadronic matrix elements of the corresponding operators eval- 
uated between a B meson and the nK isospin state with I — |. Since our intention 
here is to look for New Physics effects rather than doing precision calculations, it 
will be sufficient for our purposes to evaluate these matrix elements in a given New 
Physics scenario using the naive factorization approximation and neglecting small 
SU(3)-breaking effects. 2 Then the strong-interaction phases 4> a and 4> b vanish, and 
we obtain 

\\ u \ (a + ib) « - (cf W + ct W ) + \ ^ (cf w + ct W ) 

_ fci ^30/23 / C EW + I C EW + ^EW + I ? (3 5) 

4 V 3 3 / 

where cf w = cF w (m w ) - cf w (m w ), and 

x = = 2 < . (3.6) 

(m s + m d ) m b (m u + m d ) m b 

Parity invariance implies that the relevant hadronic matrix elements of the operators 
Oi have the opposite sign compared with those of the operators 0«. The fact that cf w 
and cf w enter with the coefficient —1 in (3.5) is a model-independent result free of 
hadronic uncertainties, irrespective of whether these coefficients receive New Physics 
contributions or not. It follows because the isovector components of the penguin 
operators O3 and Of can be related to the usual current-current operators by a Fierz 
transformation [29, 32]. For the numerical analysis we choose the renormalization 
scale /i — m b and take k = a s {m^) / a s (mw) = 1-83 and x — 1-18, yielding 

|A U | (a + ib) w -(cf w + cT) + 0.26(cf w + c™) - 1.40(cf w + c* w ) ■ (3.7) 

Next we consider the New Physics contributions to the parameter p. In the 
factorization approximation, we find 



3e 



3/2 



4|A U | 



Im 



[4 - ct -x4) + \ (4 -4-Xcf)- as 9 C 8g \ , (3. 



where cf = cf(rrib) — cf(rrib). Because the QCD evolution of the Wilson coefficients 
cf and cf is complicated, we prefer to present the result in terms of coefficients 



2 According to [24], naive factorization gives the leading term in the heavy-quark limit. 
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renormalized at the scale m b . Note that by pulling out a factor of e 3 / 2 on the right- 
hand side we avoid the difficulty of calculating the overall penguin amplitude P in 
(3.1). The contribution of the chromomagnetic dipole operator is formally of next- 
to-leading order in a s and thus could be dropped; however, we keep it because in 
some New Physics models the coefficient C% g can be enhanced with respect to its 
SM value by an order of magnitude [33, 34, 35, 36]. Using the QCD factorization 
approach of [24], we find that 

2a 

a 8g = — ^ (1 + x)(-A t ) » 4 x 1(T 3 . (3.9) 
on 

Note that the magnitude of the left-hand side in (3.8) is bounded by unity. This 
implies a nontrivial upper bound on the possible CP-violating New Physics contri- 
butions to the penguin operators, given by 



Im 



^-ci- x -4) + ^(c d 6 -4- X cf)-a 8g Cs 



4|AJ 

<J_»I. (3.10) 



Using \V ub /V cb \ = 0.085 ± 0.015, corresponding to |A„| = (7.5 ± 1.3) x 10~ 4 , and 
taking for e 3 / 2 the value in (3.3), we find that the right-hand side of this bound is 
less than 7.4 x 10~ 3 at 90% confidence level. For comparison, we note that in the 
SM the magnitude of the combination of penguin coefficients entering above is about 
3 x 10~ 3 . However, due to the smallness of the weak phase of A t the imaginary part 
of this combination is much smaller. 

3.2 Model-independent bounds on in the presence of New Physics 

The most important observable in the exploration of New Physics effects in B — > tcK 
decays is the ratio of the CP-averaged branching ratios for the decays B ± — > ^K 
and B ± — > n°K ± , given by 

where the quoted experimental value is derived from data reported by the CLEO 
Collaboration [25] . It will often be convenient to consider a related quantity defined 

as 



'It 1 ~ 1 

X R = -SL_ = 0.72 ± 0.98 cxp ± 0.03 th . (3.12) 

^3/2 

Because of the theoretical factor i?su(3) entering the definition of e 3 / 2 in (3.2) this is, 
strictly speaking, not an observable. However, the irreducible theoretical uncertainty 
in Xr is so much smaller than the present experimental error that it is justified to 
treat this quantity as an observable. The advantage of presenting our results in terms 
of Xr rather than is that the leading dependence on e 3 / 2 cancels out (see below). 
Also, some experimental errors cancel in the ratio in (3.12) [14]. 
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When writing theoretical expressions for the quantities and X R we eliminate 
the two parameters e 3 / 2 and p in favor of the measurable parameter e 3 / 2 and a "weak 
phase" y? G [-90°, 90°] defined by 

VI + p 2 V 1 + p 2 

The most direct way of probing this phase is via the direct CP asymmetry in the 
decays B ± — > 7i ± K°, which is given by 

In the SM, sin2y9 2e a sin7 is of order a few percent, and with realistic values 
for 4> p of order 10°-20° [24] one expects a very small CP asymmetry. However, in 
New Physics scenarios with new CP- violating couplings sin2<p may be significantly 
larger than in the SM. An experimental finding that A C p(ir + K°) = 0(10%) would 
constitute strong evidence for the existence of such an effect. 
The exact theoretical expression for i?* is 

i?" 1 = 1 + 2^3/2 COS if a cos(0 + 4> a ) — cos 7 COS </>J 

— 2e 3 / 2 sin ip b cos(0 + (fib — <fi P ) — sin 7 cos(0 — (fi p ) 
+ £3/2 (l — 2a cos 7 cos (fi a — 26 sin 7 cos 06 + a 2 + o 2 ) . (3.15) 

In the SM b = 0, a = o~ew, (fia = 0, and (p to very good approximation. Therefore 

iC 1 = l + 2e 3/2 (5 E w - cos 7) cos + ^3/2(1 - 25 E w COS7 + 5ew) 

< (l + 53/2 |o~ew - cos 7|) +£3 /2 sin 2 7. (3.16) 

In the second step we have used the fact that | cos 0| < 1 to obtain an upper bound on 
R* 1 . Similarly, a lower bound can be derived, which is obtained by changing the sign 
of £3/2 in the above inequality. These bounds imply nontrivial constraints on cos 7 
provided that i?* differs from 1 by a significant amount. In Figure 1, we show the 
resulting lower and upper bounds on the quantity X R versus 7, obtained by scanning 
the input parameters in the intervals 0.15 < e 3 / 2 < 0.27 and 0.49 < <5ew < 0.79. The 
latter value is obtained by using (Kfe/Kbl = 0.085 ± 0.015 in (3.4). The dependence 
on the variation of £3/2 is so small that it would almost be invisible on the scale 
of the plot. Note that the extremal values R* can take in the SM are such that 
\Xr\ < (1 + <5ew) irrespective of the value of 7. A value exceeding this limit would 
be a clear signal for New Physics [9, 17]. In view of the present large error on Xr, 
this is still a realistic possibility. Because the upper and lower bounds on X R are, to 
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a very good approximation, symmetric 
around Xr = 0, we will from now on 
only show the upper bounds, since Xr > 
is the region favored by experiment. 
Also, since the bounds change little under 
variation of <5ew and £3/2, we will work 
with the central values <5ew = 0.64 and 
£ 3/ / 2 = 0.21. By the time the experimen- 
tal data will be sufficiently precise to per- 
form the New Physics searches proposed 
in this work, the errors on these parame- 
ters are likely to be reduced by a signifi- 
cant amount. 

Let us first discuss the case where 
New Physics induces arbitrary CP-viola- 
ting contributions to the B — > nK de- 
cay amplitudes, while preserving isospin 
symmetry. Then the only change with re- 
spect to the SM would be that the weak 
obtain 




75 100 125 150 175 
gamma 

Figure 1: Upper and lower bounds on 
Xr versus I7I (in degrees) in the SM. The 
dashed line shows the bounds obtained us- 
ing the central values of the parameters e 3/ / 2 
and Sew- The band shows the current ex- 
perimental value of Xr with its la variation. 



phase (p may no longer be negligible. We 



r: 1 



2e 



3/2 



cos p (5ew — cos 7) cos (p + sin p sin 7 cos(0 — 



+ £3/2(1 - 25 EW cos7 + 5i w ) 



'3/2 

< 1 + 2e. 



3/2 



costp \8ew — COS7I + I siny9sin7| 



+ £3/2(1 - 25 E w cos7 + 5l w ) 
< ( 1 + £3/2-\/l - 2<5 E w cos 7 + <5ew 



(3.17) 



In deriving the upper bounds we have varied the strong-interaction phases <fi and 
(f) p independently. Analogous lower bounds are obtained as previously by changing 
the sign of £3/2. The last inequality in (3.17) is remarkable in that it holds for 
arbitrary isospin-conserving New Physics effects no matter how large they are. The 
corresponding bound on \Xr\ is 



Xr\ < sj 1 - 25 E w cos 7 + S EW < 1 + Sew ■ 



(3.18) 



Note that the extremal value is the same as in the SM, i.e., isospin-conserving New 
Physics effects cannot lead to a value of \Xr\ exceeding 1 + <5ew- in the left-hand 
plot in Figure 2 we show the upper bound on Xr versus 7 in the SM, and for New 
Physics scenarios with different values of ip. The three choices of <p shown correspond 
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to |p| ~ 0.27, 0.58 and 1. The gray curve shows the upper bound obtained by varying 
if. We observe that isospin-conserving New Physics can enhance the value of Xr 
relative to the SM, but only by a moderate amount. 




25 50 75 100 125 150 175 25 50 75 100 125 150 175 



gamma gamma 



Figure 2: Upper bound on Xr versus I7I for the SM (solid) and with different New 
Physics contributions. Left: Isospin-conserving New Physics with \<f\ = 15° (dashed), 30° 
(short-dashed), and 45° (dashed-dotted). The gray curve shows the bound (3.18) valid 
for arbitrary ip. Right: Isospin- violating but CP-conserving New Physics with |o| = 0.5 
(dashed), 1 (short-dashed), and 1.5 (dashed-dotted). 



Next we consider New Physics effects that violate isospin symmetry, but we first 
restrict ourselves to the important subclass of models which do not contain significant 
new CP-violating phases. Then ip and b still vanish, and we obtain 



R 



-1 



1 + 2e 3/2 



< 



a cos(0 + (pa) — cos 7 cos 

2 



+ £3/2(1 — 2a cos (p a COS7 + a 1 



1+^3/2 OI + I COS7|) 



+ £3/2 sin2 7 ■ 



(3.19) 



Again, a lower bound can be obtained by changing the sign of £3/2. In contrast to the 
previous case, now the maximal value of \Xr\ is given by 1 + \a\ and thus can exceed 
the SM bound provided that \a\ > <5ew- This is illustrated in the right-hand plot 
in Figure 2, where we show the resulting upper bound on Xr versus 7 for different 
values of \a\. In contrast with the case of isospin-conserving New Physics, even a 
moderate enhancement of the coefficient a corresponding to a 10%-20% change in 
the decay amplitudes can lead to a significant increase of the upper bound on Xr. 

If both isospin- violating and isospin-conserving New Physics effects are present 
and involve new CP-violating phases, the analysis becomes more complicated. Still, 
it is possible to derive from (3.15) a series of bounds on R~ x . We find 



R* 1 < 1 + 2^3/2 cosip (|a| + I C0S7I) + I sin<y2| + | SU17I) 
+ £3%[(H + I cos 7 |) 2 + (|6| + I sin 7 |) 2 



13 



< 



1 + e z / 2 ^(\a\ + | COS7I) 2 + (j&j + I sin7|) 1 



< 



l + e 3 /2 (l + Va 2 + 6 2 )] , (3.20) 

where in the second and third steps we have eliminated ip and 7, respectively. As 
before a series of lower bounds is obtained by changing the sign of e 3 / 2 - The corre- 
sponding bounds on Xr are 

2 



X R \ < J(\a\ + I COS7I) 2 + (|6| + I sin 7 |) 2 < 1 + Va 2 + b 2 <— + X 

£3/2 



R 



(3.21) 



where the last inequality is relevant only in cases where \/a 2 + b 2 3> 5ew- With the 
current values for e 3 / 2 and Xr, the right-hand side is less than 15 at 90% confidence 
level. The important point to note is that in the most general case, where b and p 
are nonzero, the maximal value Xr can take is no longer restricted to occur at the 
endpoints 7 = 0° or 180°, which are disfavored by the global analysis of the unitarity 
triangle [37]. Rather, the maximal value Xr 3 * = 1 + V a 2 + b 2 now occurs at 

|tan 7 | = \p\= - . (3.22) 
ct 

The situation is illustrated in Figure 3, 
where we show the upper bound on Xr 
for arbitrary New Physics contributions 
satisfying a 2 +b 2 = 1, but for different val- 
ues of a and b chosen such that the max- 
imum occurs at the endpoints I7I = 0° or 
180° (6 = 0), at the intermediate points 
I7I = 45° or 135° (a = b), and in the 
center at I7I = 90° (a = 0). The cor- 
responding values of the New Physics pa- 
rameter p required to reach the maximum 
are \p\ = 0, 1, and 00 (i.e., \p\ ^> 1), re- 
spectively. 

The present experimental value of Xr 
in (3.12) has too large an error to deter- 
mine whether there is any deviation from 
the SM prediction. If Xr turns out to be 
larger than 1 (i.e., only about one third 
of a standard deviation above its current central value), then an interpretation of 
this result in the SM would require a large value I7I > 96° (see Figure 1), which 
may be difficult to accommodate. This may be taken as evidence for New Physics. 
If Xr > 1.3, one could go a step further and conclude that this New Physics must 
necessarily violate isospin. 
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Figure 3: Upper bound on Xr versus I7I 
for arbitrary New Physics contributions sat- 
isfying a 2 + b 2 = 1. The curves correspond 
to a = 1 and 6=0 (dashed), a = b = l/\/2 
(dashed-dotted), and a = and 6 = 1 
(short-dashed). The SM bound is shown by 
the solid line. 
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3.3 New Physics effects on the determination of 7 

A value of the observable i?* which violates the SM bound (3.16) would be an exciting 
hint for new isospin-violating penguin contributions from New Physics. With the 
current central value of R* derived from CLEO data this is still a realistic possibility. 
However, even if a more precise measurement will give a value that is consistent with 
the SM bound, B ± — > irK decays still provide an excellent testing ground for physics 
beyond the SM. In the SM, the weak phase 7, along with the strong-interaction 
phase 0, can be determined up to discrete ambiguities by combining measurements 
of and an asymmetry A, which is defined as a linear combination of the direct CP 
asymmetries in the two — > ttK decay channels [8, 9]. The discrete ambiguities 
can be resolved using information on the strong-interaction phase from theoretical 
approaches such as the QCD factorization theorem derived in [24]. The theoretical 
uncertainty on the value of 7 is typically of order 10°. Although New Physics may 
not be exotic enough to lead to a violation of the general bounds derived in the 
previous section, it may still cause a significant shift in the extracted value of 7. 
This may lead to inconsistencies when the value j^k extracted in B ± — > i\K decays 
is compared with determinations of 7 using other information. 

A global fit of the unitarity triangle combining information from semileptonic 
B decays, B-B mixing, CP violation in the kaon system, and mixing-induced CP 
violation in B — > J/ipKs decays will provide information on 7, which in a few 
years will determine its value within a rather narrow range [7, 37]. Such an indirect 
determination could be complemented by direct measurements of 7 using, e.g., B — > 
DK(*> decays [38], or using the triangle relation 7 = 180° — a — (3 combined with 
a measurement of a in B — > nn or B — > np decays [7]. In our discussion below we 
will assume that a discrepancy between the "true" 7 = arg(V^* b ) and the value 7^ 
extracted in B^ — > nK decays of more than 25° will be observable after a few years 
of operation at the B factories. This will set the benchmark for sensitivity to New 
Physics effects. 

In order to illustrate how big an effect New Physics could have on the value 
of 7 extracted from B ± — > nK decays, we assume for simplicity that the strong- 
interaction phase is small. We expect that this is indeed a good assumption, since 
the QCD factorization theorem of [24] predicts that = 0(a s , A/m b ). In this case, 
a measurement of the asymmetry A will provide information about (and serve as 
a test of our assumption), whereas 7 is determined by R* alone. In the context of 
the SM, the solution obtained with cos0 ~ 1 is 



Given the present uncertainties in e 3 / 2 and <5ew> this result is a good approximation 
to the exact solution as long as |0| < 25°. Let us now investigate how New Physics 



cos ^ nK <5 EW - 




2^3/2 



(3.23) 
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Figure 4: Contours of constant Xr versus 7 for different New Physics parameters, 
assuming 7 > 0. In each plot one of the three parameters (p, a and b is varied, while the 
other two are kept fixed at their SM values. The horizontal lines indicate the SM results. 



may affect the results of this extraction, neglecting, for the purpose of illustration, all 
strong-interaction phases. As in the previous section, we focus first on the situation 
where the New Physics conserves isospin symmetry. Inserting for i?" 1 in (3.23) the 
expression given in (3.17), we obtain 

cos j %K = cos(7 + if) + <5 E w (1 - cos tp) + 0(e 3 / 2 ) • (3.24) 

For small values of if, the result is simply given by 7^ ^ 7 + y?. Therefore, to 
have a significant shift requires having \ip\ > 25°, which corresponds to rather large 
values \p\ > 0.5 and hence an 0(1) change in the decay amplitudes. The situation 
is different if the New Physics contributions violate isospin symmetry. Consider, for 
simplicity, the case where there are no new CP- violating phases. From (3.19) it then 
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follows that 

cos 7^ = cos7-a N p + 0(£3/ 2 ) , (3.25) 

where 

a NP = a — 5 E w ■ (3.26) 

Now even a moderate New Physics contribution to the electroweak penguin coeffi- 
cients can lead to a large shift in 7. In the most general case, where all New Physics 
contributions are present, we obtain 

cos 7^ = cos(7 + <p) + 5ew - acostp + bsimp + 0(e 3 / 2 ) , (3.27) 

which again allows for large shifts. 

These observations are illustrated in Figure 4, where we show contours of con- 
stant Xr versus 7 for different values of one of the parameters ip, a and b, keeping the 
other two fixed to their SM values. Without loss of generality we assume that 7 > 0. 
These plots show that even moderate New Physics contributions to the parameter 
a can induce large shifts in 7. On the other hand, small values of ip lead to much 
smaller effects. Likewise, the effects induced by the parameter b are much smaller. 
We stress that the present central value of Xr « 0.7 is such that negative values of 
ip, as well as values of a less than the SM result a ~ 0.64, are disfavored since they 
would require values of 7 exceeding 100°, in conflict with the global analysis of the 
unitarity triangle [7, 37]. 
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Figure 5: Difference j n K — 7 versus Xr for different New Physics scenarios, assuming 
7 > 0. Left: a = <5 E w, b = and <p = 0, ±15°, ±30°, and ±45°. Right: <p = b = and 
inp/Sew = 0) ±0.5, ±1, and ±1.5. Solid curves refer to positive values of <p and onp, 
dashed curves to negative ones. 

In Figure 5, we show the difference 7^ — 7 as a function of Xr for the two 
cases where either ip or a are varied with respect to their SM values. The dashed 
curves in the plots refer to negative values of tp or a^p, which are disfavored by the 
present value of Xr. This implies that isospin-conserving New Physics can only lead 
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to moderate shifts in 7, which reach the 30° level for large values p = 0(1). Isospin- 
violating New Physics effects, on the other hand, can induce very large shifts of 7 even 
if they are of moderate size. As an example, consider the case where a future, precise 
measurement would yield Xr = 1. An interpretation of this result in the SM would 
imply a relatively large value 7 ~ 106°, which can be determined with a theoretical 
uncertainty of about 10° or better [8, 9]. Imagine that all other information about the 
unitarity triangle favors a value of 7 ~ 75°, again with a small error. To accommodate 
this difference, one could either invoke a large, isospin-conserving but CP-violating 
New Physics contribution such that ip ~ 45° (corresponding to p ~ 1), or an isospin- 
violating electroweak penguin contribution such that a is twice as large as in the SM. 
The first solution would imply a New Physics contribution to the decay amplitudes 
of order 100%, whereas the second would imply only a small contribution of less than 
15%. 

3.4 "Wrong kaon" decays 

So far, when considering B — > 11K decays we have implicitly assumed an underlying 
quark transition of the form b — > sqq, in which case the decays with a neutral kaon 
in the final state are B + — > ir + K° and B~ — > ir~K°. Indeed, in the SM this is an 
excellent approximation, because the quark transition b — > dsd leading to the "wrong 
kaon" decays B + — > ir + K° and B~ — > ir~K° is highly suppressed. However, this may 
no longer be the case in the presence of New Physics. 3 

In practice, only the B^ — > tt^Ks^l decay rates can be measured. In particular, 
the CLEO result for quoted in (3.11) really refers to the ratio 

^* " B (5+ - tt0^+) + B(5- - vro^-) " + ' rong ' j ' [6 - M) 
which differs from by the "wrong kaon" contribution 

\Y | 2 _ B(g + ^7r+^) + B(g-^vr-^) 

' Wr0nS ' B(fi+ ^ TT + ^0) + B(S" ^ TT-^0) • 1 J 

Note that the presence of this contribution could only enhance the observed value 
Rl xp with respect to i?*. This observation, combined with the fact that R cxp is not 
much larger than the expected value for R* in the SM, allows us to put bounds on 
the "wrong kaon" contribution in specific New Physics models. 

In analogy with (2.5), there are three operators entering the effective Hamiltonian 
for b — > dsd decays, which we define as 

Oi d = (b a d a )v-A(si3dp)v+A , 
= (b a dp)v-A(spd a )v+A , 
Of = (b a d a ) v _ A {spd p ) v „ A , (3.30) 



3 For the related decay b — > sds, this possibility has been explored in [39]. 



18 



and similarly there may be operators dd of opposite chirality. We denote the corre- 
sponding Wilson coefficients by c dd and c dd , respectively. The renormalization-group 
evolution of these coefficients can be read off from (2.11) and (2.12) by replacing 
c ew _^ c dd^ c ew _^ c dd ^ an j c ew _^ c dd \j s [ n g factorization, and normalizing the 

result to £3/2, we find 



Y - -^1 

^ wrong I > I 



-dd 2 3/23 ^d , %X 1 ^30/23 {-dd , 1 -dd 
c 3 3 K 1 ^ K \ c 2 "r g t-i 



(3.31) 



where cf d = c dd {mw) — c dd {m w ). Inserting this result into (3.28), and performing 
the evolution to the scale /i — m&, we obtain 



cf - 0.26cf + 1.40^1 « piiW^- - 1 . (3.32) 




Using the current values of A M , £ 3 / 2 and i?° xp , and taking > 0.45 corresponding 
to the smallest possible value in the SM, we find that the right-hand side is less than 
5.2 x 1(T 3 at 90% confidence level. 



4. Trojan penguins from tree-level processes 

In this and the following sections, we calculate the matching conditions for the pen- 
guin coefficients cf and cf in a variety of extensions of the SM with new flavor- 
violating couplings. We focus first on models where such couplings appear at tree 
level. Loop- mediated New Physics contributions will be discussed in Section 5. 

4.1 Flavor-changing Z-boson exchange 

A generic feature of extensions of the SM with extra nonsequential quarks is the 
presence of tree-level flavor- changing couplings of the Z boson, such as a bsZ vertex. 
A detailed discussion of such models can be found, e.g., in [40]. The relevant terms 
in the Lagrangian read 

£fcnc = -t-^tt E * [«2V(l - 75) + + 75)] d% , (4.1) 

4costV i _ L . 

where i,j are generation indices. The quantities K b £ and K b £ are the two new complex 
parameters relevant to b — > s transitions. Since the flavor-violating interactions are 
small (see below), the flavor-diagonal couplings of the Z-boson are to leading order 
the same as in the SM. It follows that at low energies the tree- level Z exchange for 
the decay b — > sqq leads to the effective Hamiltonian 

H cS = % K s (h) v _ A + k%(B8)v + a] E [Cl(qq)v-A + C%(qq) v+A ] , (4.2) 
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where 



a = l -- 2 -sm 2 e W) c u R '^hro 



C t = -\ + \^ 2 e w , C d R = ^sm 2 9 w . (4.3) 

It is straightforward to match this result with the generic form of the penguin terms 
in the effective Hamiltonian, and to deduce the corresponding values of the Wilson 
coefficients. The nonvanishing coefficients at the matching scale are c\ = K b £ C R , 
c| = n b £ C q L , c? = 4 C q L , and 2| = n b £ C q R . In the notation of (2.7), this implies 

c QCD = _^f ? g QCD = _^ ) c QCD = g QCD = 0; ^ 

as well as 

C ? W = -4W0 W , C f W = K b L S COS 2 6 W , 

cf w = 4 s cos 2 ^vk , 2f W = -4 sin 2 ^ . (4.5) 

These electroweak penguin coefficients would be of the same order as the SM result 
for cf w in (2.10) if |4 s .rI ~ 5 x 10~ 4 which, as we will see below, is consistent with 
experimental bounds. Inserting these results into (3.5) and using |A„| ~ 7.5 x 10~ 4 
yields 

a N p + ib « -1.1 x 10 3 (4 + 0.51/$) • (4.6) 

We stress that the simple model considered here is a prototype of New Physics models 
in which the electroweak penguin coefficients cf w and cf w are not suppressed relative 
to the QCD penguin coefficients c^ CD and c^ CD . This property is in contrast with 
the SM, where electroweak penguins are suppressed by small gauge couplings. 

At present, the strongest constraints on n b £ and 4 follow from the experimental 
bound on the B — > X s e + e~ decay rate. Since this bound lies far above the SM 
prediction for this process, we can neglect the SM contribution and write the effective 
Hamiltonian for this process as 

H eS = ^ [K»(b8)v-A + *%(b8)v + A] VI {ee) V -A + C R (ee) v+A ] , (4.7) 

where 

C£ = -- + sm 2 9 w , C e R = sm 2 9 w . (4.8) 

It is convenient to normalize the result for the B — > X s e + e~ decay rate to the 
semileptonic rate. Then many common factors cancel, and we obtain 



T(B ^ X s e + e-) 


K bs 


2 +l 


K R\ 


2 


T(B -> X c e~v e ) f(m c /m b ) \ V cl 


|2 
| 



(Ct) 2 + (C R ) 2 } ^157(|4| 2 + l4l 2 ) , (4-9) 
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where we have used jV^I 0.04, and f(m c /m b ) m 0.5 for the phase-space factor 
in the semileptonic decay. Using the upper bound B(B — > X s e + e~) < 5.7 x 10~ 5 
together with B(B -> X c e-z/ e ) = 0.105 [41] yields 

\I\ k l\ 2 + \ k r\ 2 < 19 x 10 ~ 3 • ( 4 - 10 ) 
Combining this result with (4.6), we obtain the upper bound 

y/a$, p + b 2 < 2.0 , (4.11) 

which may be compared with the SM value a ~ 0.64. It follows that tree-level Z 
exchange with new flavor-violating couplings can yield isospin-violating electroweak 
penguin effects that are up to a factor 3 larger than in the SM. 

For completeness, we also mention the resulting bound on the New Physics 
parameter p. Neglecting renormalization-group effects, we obtain from (3.8) the 
estimate 



-'29 



0.8Im4 s + Im/4 S =>■ |p|<0.05, \<p\ < 3° , (4.12) 



VT+7 

indicating that in this model there is no room for large values of p. 
4.2 Extended gauge models with a Z' boson 

A new neutral boson Z' with tree-level flavor- changing couplings to quarks is a generic 
property of many models with an extended gauge group. The analysis of electroweak 
penguins in such models is very similar to the flavor- changing tree-level Z exchange 
discussed above. For simplicity, we assume no significant mixing between the Z and 
Z' bosons. Then the effective Hamiltonian is a simple generalization of (4.2), i.e. 

Waff = 9 -^fWt(bs)v-A + 4 S (bs) v+A ] £ [C q L (qq)v-A + C q R {qq) v+A } , (4.13) 

m Z' q 

where C\ R now denote the charges of the quarks under the new U(l)' group. Intro- 
ducing the ratio 

and taking into account that, since we neglect Z—Z' mixing, the Z 1 charges are the 
same for all left-handed fields (in particular, C£ = Cf), we find 

QCD _ , /bs U R + lL, R „QCD _ tJbsr* 

c? CD = i^Cl , c? CD = £4* ^±^i , (4.15) 
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and 



cf w = i4*{C u R - C d R ) , c™ = £K% s (C«-C d R ), cf w = cf w = 0. (4.16) 
Inserting these results into (3.5) gives 

a NP + ib » ^ C R~ C *) (4^ + 0.26/cf) . (4.17) 

The allowed range for the relevant parameters in Z' extensions of the SM is 
largely model dependent. For example, the bounds derived from the upper limit on 
the B — > X s l + l~ branching ratio depend on the lepton charges under the new U(l)' 
gauge group. In the so-called "leptophobic" Z' models these charges are arranged so 
as to vanish or be very small [42]. Therefore, in general the contributions of the Z' 
couplings to the electroweak penguin operators can be arbitrarily large. In fact, the 
best model-independent bound on these couplings follows from the second inequality 
in (3.21), which implies 

\aC u R - C d R )\ |4 S + 0.264 s | < 0.01 , (4.18) 

where we have neglected the small SM contribution. Furthermore, assuming C q LR = 
0(1) and no cancellations, the bound (3.10) gives |£Im4 S R| < O(10~ 2 ). Turning 
these observations around, we conclude that in extended gauge models with flavor- 
changing Z' couplings such that £4fi? = O(10~ 2 ) there can be very large New Physics 
effects in B — > ixK decays. 



4.3 SUSY models with R-parity violation 

In SUSY models with broken R-parity extra trilinear terms are allowed in the su- 
perpotential, some of which can give rise to a large enhancement of the electroweak 
penguin coefficients. Denoting by L % L) Q l L , u R and d l R the chiral superfields contain- 
ing, respectively, the left-handed lepton and quark doublets, and the right-handed 
up- and down-type quark singlets of the i-th generation, these terms read 

W = K jk Ll Q j L 4 + Kjk u R 4 4 • (4-19) 

At low energies, slepton and squark exchange can generate local penguin operators. 
The most general case has been treated in [18]. For simplicity, we neglect left-right 
sfermion mixing, which is a small effect and does not generate new operators. We 
then find for the coefficients of the various penguin operators 

3 \/* \/ \//* \ll 

~u _ A »12 A il3 ~u _ pi _ A H3 A 112 (a 9f| N 

° 2 ~ h ^G F ml L ' C3 " C4 - 2V2G F m^ ' ^ 
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and 



3 \'*V 3 \'*V 

d _ ST A »31 A t21 d _ A i32 A ill 



3 \/* \/ 3 V* V 

~d _ A il2 A il3 ~d _ A «ll A i23 



_^ = g^ = V C! 13Aa2 9 . (4.21) 

Previous authors have investigated bounds on some of these R-parity violating cou- 
plings in the context of nonleptonic B decays [43, 44]. However, in these studies 
model-dependent predictions for the overall penguin amplitude P in (3.1) are em- 
ployed. The only significant bound which has an impact on our analysis comes from 
a combination of constraints derived from limits on double nucleon decay into two 
kaons and neutron-antineutron oscillations, yielding lA'/j^A'/^l < 10~ 9 [43]. There- 
fore, is it safe to neglect c 3 and c\. 

Not all of the above coefficients contribute to the isospin- violating terms parame- 
trized by a and b. Up to a small SU(2)x, breaking in the slepton and sneutrino masses 
we find = c 2 , and thus cf w ~ 0. Moreover, in (3.5) only the sum cf w + cf w 
contributes, which vanishes since c| = — c|. It follows that 

a NP + ib « 2.83 x 10 3 £ (10 ° G 2 eV)2 (Af 31 A; 21 + A^Ajn - A^A^) . (4.22) 

i=l m vi 

The result for the parameter p is more complicated. Setting for simplicity all sfermion 
masses equal, we find from (3.8) 

P ~i06 (10 ° GeV)2 



Vi + p 2 m ) 

x £ imUA^A:; 2 + 3 X (a:i 2 a: 13 - a: 31 a: 21 ) + - . (4.23) 



Using the results derived in Section 3, we can obtain bounds on several of the R-parity 
violating couplings. Assuming that only one combination of couplings is dominant, 
neglecting the SM contribution, and using a common sfermion reference mass of 
100 GeV, we find from (3.21) that at 90% confidence level 

IE A5iA; 21 <4.9x10" 3 , 



i=i 

3 



£Af 32 A; n | <4.9xl0- 3 , 
i=i 

<4.9xl0- 3 . (4.24) 
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In addition, from (3.10) we obtain 



£lrn(A^ 2 )| <3.7xKT 3 , 
EM^i^ai) <4.1xl(T 3 , 



i=i 

3 



i=l 



£lm(A£ 2 Ay <4.1xl0" 3 , (4.25) 



where we do not present bounds on the imaginary parts of couplings which are weaker 
than the corresponding bounds on the absolute values in (4.24). 

It is interesting that SUSY models with R-parity violation provide an example 
of scenarios in which b — > dsd transitions may not be suppressed relative to b — > sdd 
transitions. As pointed out in Section 3, this can lead to potentially large "wrong 
kaon" decays of the type B + — > tt + K° and B~ — > tt~K°. In the model considered 
here, the only nonvanishing coefficients are 

3 w* w 3 \/* \i 

c f = V 7 lAtl2 9 , cf = V 7 lAa3 9 . (4.26) 

The result (3.32) can be used to obtain the bounds 

lASiAJul < 3.4 x 10- 3 , |A^Ay < 3.4 x 10~ 3 , (4.27) 

again at 90% confidence level. 

Our bounds in (4.24) and (4.25) are stronger than the ones discussed in the 
literature [43, 44] and refer to a larger number of i?-parity violating couplings. Most 
importantly, however, they are affected by much smaller hadronic uncertainties. The 
bounds in (4.27), on the other hand, are weaker than constraints derived from B-B 
and K-K mixing [45]. 



5. Trojan penguins from loop processes 

Having considered in the previous section some specific models with tree-level FCNC 
couplings, we now explore extensions of the SM in which new contributions to the 
electroweak and QCD penguin operators arise at one-loop order. In particular, we 
study in detail the structure of electroweak penguins in SUSY models, where isospin- 
violating b — > sqq transitions can arise due to strong-interaction gluino box diagrams. 
This provides another realization of models in which electroweak penguins are not 
suppressed relative to QCD penguins. For completeness, we also discuss two-Higgs- 
doublet models and models with anomalous gauge-boson couplings. They are simple 
since there are no new CP-violating phases, so only the parameter a can receive a 
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New Physics contribution. However, there is no parametrical enhancement of the 
electroweak penguins relative to the SM, and thus the New Physics contributions 
tend to be small. 



5.1 SUSY models 

In SUSY extensions of the SM with conserved R-parity, the potentially most impor- 
tant contributions to the Wilson coefficients of the penguin operators in the effec- 
tive Hamiltonian (2.6) arise from strong-interaction penguin and box diagrams with 
gluino-squark loops. They can contribute to FCNC processes because the gluinos 
have flavor-changing couplings to the quark and squark mass eigenstates. Provided 
there is a significant mass splitting between the right-handed up and down squarks, 
gluino box diagrams are also the most important source of isospin violation in SUSY 
models [15]. The corresponding contributions to the electroweak penguin coefficients 
are then much more important than other SUSY contributions from photon or Z 
penguins usually discussed in the literature. In fact, in such a scenario SUSY contri- 
butions to the coefficients of the electroweak and QCD penguin operators are of the 
same order and scale like a"j. / m^ VSY , where msusY is a generic mass of the superpar- 
ticles. Whereas the QCD penguin contributions are typically smaller than in the SM, 
the electroweak penguin contributions can be important, since their scaling relative 
to the SM coefficients is controlled by the ratio (a s /a)(m^/mg USY ) ~ 1. In our 
analysis we will consider only these potentially large gluino box and penguin contri- 
butions and neglect a multitude of other SUSY diagrams, which are parametrically 
suppressed by small electroweak gauge couplings. The latter include photon or Z 
penguins with gluino-squark, chargino-squark, neutralino-squark, or charged-Higgs- 
quark loops, and various box diagrams containing at least one chargino or neutralino. 
We have calculated all of these diagrams and, for generic regions in SUSY parameter 
space, have found their contributions to be largely suppressed relative to the pure 
gluino diagrams. 

Large SUSY contributions to the penguin operators via gluino loops require near 
maximal mixing between the strange and bottom squarks, so that the squark mass- 
insertion approximation is not valid. We therefore present our results using the 
general vertex-mixing method, summing over diagrams with different squark mass- 
eigenstates in the loops [46, 47]. We denote by T Dl the rotation matrices relating 
the left-handed down-squark interaction states in the quark mass-eigenbasis, q l L (q = 
d, s, b), to the squark mass eigenstates, di (i = 1, . . . , 6), such that q\ = (T® L )* d iy 
with obvious generalizations for the up- and right-handed squarks. In addition, 
Xg t g = ra~./m~, where is the mass of the i-th down (q = d) or up (q = u) 
squark mass eigenstate. In the operator basis of (2.5), we obtain for the gluino box 
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^6, box 



contributions to the Wilson coefficients at the SUSY matching scale 

c l,box = 2^G F m 2 ~ ^ ( r i/) r jcf ^8 F ( X d i g^ X d j g) ~ lg G ( X d i g^ X d j g) 

c 2,box = i^plGpm 2 - r ^ ( r i«*) r i«* 5 F ( x di^ x ujg) + g G i x di^ Xa jg ) 

c 2,box = 2^/2G F m 2 ~ ( r .j<f ) 5 F ^ x d^ X d j9 ) + g G ( X dig' X d j g) 

c 3,box = 2^G F m 2 - ' r * s V r i"J _ g F ( x dig' x ujg) + 3^ ^(^g' 

C 3,box = 2^j2G F m 2 ~ { ^ ( r ^) ^ _ 9 F ^ X d i g^ X d j g) + 3g G ( X d i g^ X d j g) 

+ ( r i6 L ) ^ ( r id L ) 3 F ( X d i g^ X d j g) ~ G ( X dig' X djg) 

c 4,box = 2\f2G F m 2 - r ^ ( r i« ) r i« 3 ^(^(iifl'^s) _ ^ G ^ x d\~gi x ^i9i 

Ca ^~ 2V2G F m 2 g \ [ilb ) L " ^ jd ) L > d 

+ (r^)*rf/(rf/)*rf/ 

C 5,box = 2^2G F m 2 - ^i^) ^ 18 F ^ X d i g^ X d j g) - Jg" G ( X d l g^ X d j g) 

2V2G F m 2 ( r * L ) F ? ^i^) ^ 6 F ( X d i ~g^ X d j g) + g ^(^g' 



1 7 

3 F ( X d ig i X d j g) ~ J2~ G ( X d t gi X d J g) 

5 f 

~g F ( X d t g' X djg) + gg G \ X digi X djg) 



' X dfg) 



, (5.1) 



where repeated indices are summed over, and C5 box = Cg box = 0. The functions 
F(x,y) and G(x,y) are given by [46, 47] 

a; In a; y In y 1 



F(x,y) 
G(x,y) 



(x - y)(x - l) 2 (y-:r)(y-l) 2 (x-l)(y-l)' 



x 2 Inx 



y 2 my 



+ 



1 



(5.2) 



(x - y)(x - l) 2 (y - x)(y - l) 2 (x-l)(y-l)' 

The corresponding expressions for the coefficients cf of the opposite-chirality oper- 
ators are obtained via the exchange L <-> i? in the expressions for cf. In practice, 
c 5,box; c 6,box an d ^e second terms in cf box , C4 box (as well as the corresponding terms 
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in the coefficients of the opposite-chirality operators) can be neglected due to B d -B d 
and K-K mixing constraints on the off-diagonal (1-3) and (1-2) entries of the T Dl 
and Y Dh matrices. Gluino box diagrams which in the mass-insertion approximation 
would contain left-right squark-mass insertions can also be neglected and have not 
been included in (5.1). Specifically, these graphs would contain the mass insertions 
5m 2 r , 5m 2 t , 5m 2 3 t or 5m 2 7 t , whose magnitudes are tightly constrained by the 

sl^r' s R b L ' d L b R d R b L ' & J 

experimental value of the B — > X s d 7 branching ratio [48]. Further suppression of 
such graphs can be expected on theoretical grounds since the remaining left-right 
squark-mass insertions they would contain are suppressed by light quark masses in 
general supergravity theories [49, 50] and in SUSY theories of flavor [51, 52, 53]. 

In addition to the box contributions, the QCD penguin coefficients also receive 
contributions from gluon penguin diagrams containing gluino-squark loops. These 
are given at the SUSY matching scale by [54, 55] 



r q 

l,peng 



^3,peng 



^2,peng 



pcng 



" 2V2G F m 2 ^ 12 A{X ^ ] + 9 B ^ 



(5.3) 



where 



1 (1 + 2x) hax 

{X) ~ 2(1 -x) + 6(1 -x) 2 ' 

11 - 7x + 2x 2 \nx 

B(x) = ; ; — . (5.4) 

V ; 18(1 -x) 3 3(1 -x) A V ; 

The opposite-chirality contributions are again obtained via the substitution L — > R. 

FCNC constraints on the off-diagonal entries of the squark mass matrix allow 
for a simple parametrization of the gluino box and penguin contributions to b — > suu 
and b — > sdd transitions, up to small corrections which have a negligible impact 
on the resulting Wilson coefficients. Let us first consider the down-squark sector. 
Constraints from B^-B^ and K-K mixing imply that, to good approximation, the 
down squark is decoupled from the strange and bottom squarks. We also neglect the 
left-right down-squark submatrix, since even in the most general case of supergravity 
theories with arbitrary Kahler potential its entries are much smaller than the typical 
squark-mass squared [49, 50] , 4 and in SUSY theories of flavor its entries are even more 
suppressed [53]. The above simplifications essentially give three "left-handed" and 
three "right-handed" down-squark mass eigenstates, obtained by diagonalizing the 
left-left and right-right squark submatrices. Specifically, in the physical down-quark 



4 The largest entries in the absence of flavor symmetries are suppressed by a factor mf,/msuSY- 
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basis (d,L,SL,bi) the left-left submatrix takes the form 



( m\ 


v o 



o \ 



rn 



22 



K 3 ) ! 



m 23 
m 33j 



(5.5) 



Let us denote the left-handed mass eigenstates by qi (q = d, s, b) and their masses 
by rri~ L , making the identification m~ = m\ x . Then the left-handed squark mass- 



eigenstates take the form 



d,L = di 



(1\ 



w 



sl = d 2 



V 



o 

cos L 
sin 9 L e~ 



( 



br= d a 



) 





sin 9 L e lS 
cos9l 



J 



(5.6) 

where 5l is a new CP-violating phase. We take \9l\ < 45°, so that the squark mass- 
eigenstate §l is more closely aligned with the s quark, and b^ with the b quark. 
In the case of the box graphs we also need to consider the up-squark sector. D-D 
mixing bounds [52] imply that, to good approximation, the up squark is decoupled 
from the charm squark in the sense that including the phenomenologically allowed 
mixing between the two will lead to negligible modifications of the Wilson coefficients. 
Without loss of generality, we can also ignore mixing between the up and top squarks, 
which is a good approximation in SUSY theories of flavor. 5 Finally, as before we can 
neglect the mixing between the left- and right-handed up squarks. 
Taken together, the above approximations imply 



r 



1 id — U ' 



( r D L y F D L _ c _ ( r u L y r u L _ r _ 

{ L jdJ 1 jd —°d j d L ^ \ L ju ) 1 ju —OujU 



u L ■ 



(5.7) 



The diagonalization in the right-handed sector proceeds in a similar way, leading to 
mass eigenstates parameterized by a mixing angle Or and a weak phase 5r. 

It is now straightforward to reexpress the gluino box and penguin contributions 
in (5.1) and (5.3) in terms of our parametrization. The combined results for the 
coefficients c- read (q = u, d) 



a 2 s sin 20 L e i&L 
4 v / 2G F m? 
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5 12 

G(x~ bL ~ g , Xg R g) + - A{x~ bL ~) + - B{x- b J 



18 



5 Although up-top squark mixing can be large in supergravity theories with arbitrary Kahler 
potential, it would not modify our conclusions qualitatively. Furthermore, a model which would 
admit large up-top squark mixing but satisfy the D-D constraint on up charm squark mixing 
would have to be very contrived. 
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r q 



( X b L g ~ ^Sig) , 

g F ( x b L ~g, x<j R9 ) + « Xg R§ ) - - A{x~ bL ~) - - B{x b J 



4y/2G F ml 



C 4 



\ X b L g ~~ ^ X s L g) i 

o? a sin 26 L e iSL 
4 v / 2G F m| 

o£ sin 26 L e i&L 



1 



x 9lq) + 3^ + o A ( X b L g) + n B ( X b L g) 



9 * 



2 ^b L g! -r g ^V^^gy 



4\72GW 



I 7 3 2 

3 x «l~ 9 ) - ^ G(x~ bL§ , Xf L g) - - A(x b J - - B(x b J 



x b L h ~^ x s L g) ■ 



(5. 



The coefficients cf are obtained by substituting L <-> R above. The coefficients c| 6 
and c| 6 vanish in the approximation (5.7). 

Let us identify those regions of SUSY parameter space which can give large 
contributions to the Wilson coefficients cf and cf. It turns out that a small gluino 
mass is favored for all of the coefficients. Large contributions also require m~ L R ^> 
m~ and small ra~ , or vice versa. Both options, m\ m\ or m~ ^> m\ , 
are equivalent as far as the magnitudes of the new contributions to the Wilson 
coefficients are concerned. Perhaps the first option is more attractive given that 
constraints from K-K mixing are more stringent than those from B-B mixing. 
In addition, in models where SUSY is broken at high energies, e.g., at the grand- 
unified theory (GUT) or Planck scales, renormalization tends to make the third- 
generation squarks the lightest at the weak scale. Large contributions also require 
near maximal mixing between the strange and bottom squarks, i.e., |sin2# L | or 
I sin 26 r\ not far below 1. For the left-handed squarks this condition poses an obstacle 
for model-building due to the requirement that the CKM mixing-angle hierarchy 
must be reproduced. However, a large mixing between the right-handed strange 
and bottom squarks poses no such problem. In fact, several SUSY models of flavor 
utilizing horizontal symmetries exist in which this situation is realized [52, 56]. It 
therefore appears unlikely that SUSY contributions to the SM operators could be as 
large as those to the opposite-chirality operators. Nevertheless, to be fully general we 
will present numerical results for large left-handed or right-handed mixing. Finally, 
as mentioned at the outset, significant contributions to the isospin- violating operators 
require a large mass splitting between the up and down squarks of the first generation. 
This is possible only in the right-handed sector, since SU(2)x, invariance implies that 
to- = rn\ L up to tiny SU(2)i-breaking corrections. Therefore, only cf w , cf w , cf w 
and cf w can acquire significant gluino box contributions. The magnitudes of these 
contributions are symmetric under interchange of m\ and mh . One can consider 
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m~ 3> m? and small m? , or vice versa. 
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Figure 6: Left: QCD penguin coefficient |c^ CD | in units of 10~ 4 |sin2#£| versus the 
common mass rriu R = toj^ (left axis) and the mass splitting mg L — (right axis), 
for mi_ = nig = 250 GeV. Right: Electroweak penguin coefficient |cf w | in units of 



10 4 1 sin 28 l\ versus the mass splittings mi 



m. 



« R 



(left axis) and mg^ — mg (right axis), 



for ma R = m- bL 



rrir, 



250 GeV. All masses are given in GeV. 



We are now ready to present our numerical results. We begin with the SUSY 
contributions to the penguin coefficients at the SUSY matching scale, which for 
simplicity we take to be trw, thus ignoring the slow running of a s and superpartner 
masses above the weak scale. We find that the QCD coefficients obey the approximate 
scaling relation c^ CD ~ CD ~ — 3c^ CD ~ — 3cg CD (provided we use the same 
masses for left- and right-handed squarks), which according to (5.3) is exact for the 
contributions of the penguin diagrams but only approximate for the box diagrams. 
The coefficients cr are the same as the c? CD if all labels L *-> R are interchanged. In 
the left-hand plot in Figure 6 we show the largest coefficient, c^ CD , for a common mass 
mi L = rrig = 250 GeV as a function of the mass splitting between the left-handed 
strange and bottom squarks and of the common mass = fn^ R . Note that only 
the box contributions in (5.1) depend on the latter two masses, but not the penguin 
contributions in (5.3). We find that these two contributions interfere destructively. 
For small up- and down-squark masses the box contributions are dominant, whereas 
for large masses the boxes decouple and the penguin contributions dominate. For 
intermediate masses there is a region with large destructive interference, where c^ CD 
vanishes or takes small values. Note that typical values of the coefficients c^ CD are 
of order few times 10~ 4 (provided the gluino is as light as 250 GeV and there is 
sufficient mass splitting between the strange and bottom squarks) , which is an order 
of magnitude less than the typical size of QCD penguin coefficients in the SM [10]. 
If the gluino mass is increased and all mass ratios remain the same, then the SUSY 
contributions to the Wilson coefficients decrease, scaling like (250 GeV/rrig) 2 . 
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We next turn to the coefficients of the electroweak penguin operators, which only 
receive contributions from the gluino box diagrams. As mentioned above, because of 
SU(2)l symmetry only cf w , cf w , cf w and cf w are important, and we find that for 
equal strange-bottom mixing and mass splitting in the left-handed and right-handed 
squark sectors they roughly scale according to cf w ~ cf w ~ — 3cf w ~ — 3cf w . In the 
right-hand plot in Figure 6 we show the largest coefficient, cf w , as a function of the 
mass splittings between the left-handed strange and bottom squarks and between 
the right-handed up and down squarks. Provided both splittings are significant, 
the typical values of the electroweak penguin coefficients are of order few times 
10~ 4 , which is comparable with the value of the coefficient cf w in the SM, given in 
(2.10). Therefore, in certain regions of SUSY parameter space there can be important 
isospin-violating contributions to the parameters a and b entering the — > tcK 
decay amplitudes. An important point to notice is that SUSY contributions to the 
electroweak penguin coefficients are typically of same order as, or can be larger than, 
the contributions to the QCD penguin coefficients, if there is sufficient mass splitting 
between the right-handed up and down squarks. 

The SUSY contributions to the pa- 
rameters a and b can be decomposed as 

(a N p + ^)susY (5.9) 
= A L sin 26 L e i&L + A R sin 26 R e i&R , 

where receives contributions from the 
electroweak penguin coefficients cf w and 
cf w , and Ar receives contributions from 
cf w and cf . As previously mentioned, 
to obtain large values of these parame- 
ters requires a significant mass splitting 
between the right-handed up and down 
squarks, as well substantial splitting 
between the left- or right-handed strange 
and bottom squarks. Exchanging strange 
and bottom or up and down squarks leaves 
the results invariant up to a sign. In 
Figure 7 we show the values of Al and 
Ar versus mg R — for two choices of the strange-bottom splitting, such that 
m s L R/ m b L R = 2 or 4. We note that for large mass splittings the magnitude of A L 
can be up to twice the SM parameter 5ew ~ 0.64, whereas for more moderate split- 
tings \A L \ ~ 5ew can be obtained. The magnitude of the parameter A R is typically 
smaller by a factor of 3. According to Figure 5, SUSY contributions of this size could 
lead to shifts of up to ±50° in the extracted value of 7^. Even in the more realistic 
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Figure 7: Quantities Ar (two upper 
curves) and Al (two lower curves) versus 
the mass splitting — m,u R , for ras H = 
mr = m~ n = 250 GeV. The solid curves 
refer to rris L R = 1000 GeV, the dashed ones 
to mz, „ = 500 GeV. 
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case where only right-handed strange-bottom squark mixing is large, shifts of up to 
±25° are possible. 

Let us now briefly discuss SUSY con- 
tributions to the parameter p describing 
CP-violating but isospin-conserving New 
Physics effects in B ± — > ttK decays. From 
(3.8) it follows that QCD penguin coeffi- 
cients of order few times 10~ 4 can only 
lead to rather small values \p\ ^ 0.1. In 
analogy with (5.9), we define 

(siny?)susY (5.10) 
= S l sin 29 L sin 5l — Sr sin 29 r sin 5r , 

where Sl (Sr) depends on the mass split- 
ting between the left-handed (right-han- 
ded) strange and bottom squarks, and 
both quantities depend on the masses of 
the left- and right-handed down squarks. 
In Figure 8 we show the values of these 
quantities versus the common down-squark mass for two choices of the strange- 
bottom splitting. We see that, indeed, typical values of Sl and Sr are of order 0.1 or 
less. The corresponding contributions to p are of the same order for maximal weak 
phases and large strange-bottom mixing. According to Figure 5, SUSY contributions 
of this size can only lead to insignificant shifts in the extracted value of j % k- 

Thus far we have only considered contributions to p due to the four-quark pen- 
guin operators. The largest possible contributions in fact arise if the coefficient Cs g 
of the chromomagnetic dipole operator, or the coefficient Cs g of the corresponding 
dipole operator with opposite chirality, have very large magnitudes compared with 
their values in the SM, implying enhanced b —>■ sg transitions. This scenario has 
been discussed in the context of the low semileptonic branching ratio and charm 
yield in B decays [34, 35, 57] and the large B — > X a T]' branching ratio [58, 59, 60]. 
In SUSY models it is most easily realized via gluino-squark loops containing left- 
right strange-bottom squark mass insertions [33, 35, 36, 46]. Constraints on these 
graphs from B — > X s ^ decays allow for B(B — > X sg ) < 10%, which corresponds to 
(|C 8ff | 2 + iCgJ 2 ) 1 ' 2 ~ 1 (at the scale m b ), together with possibly large CP- violating 
phases in these coefficients [61]. From (3.8) it follows that large values p = 0(1) can 
be obtained in such a scenario. According to Figure 5, in this extreme case large 
shifts in the value of 7 vr ^- caused by isospin-conserving New Physics are not excluded. 

At present, there are no significant phenomenological constraints on the angles 
9l and 6r parameterizing the mixing between the strange and bottom squarks. In 
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Figure 8: Quantities Sl,r versus the 
common mass mr = mi , for mr = 

O.L OR> b LyR 

rrig = 250 GeV. The solid curve refers to 
wij t R = 1000 GeV, the dashed one to 
m SL R = 500 GeV. 
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particular, the measured B — > X s 7 branching ratio does not impose a useful con- 
strain on these parameters. However, an important constraint would emerge if in the 
future B s -B s mixing were found to be consistent with, or not much larger than, its 
predicted value in the SM. For simplicity, we assume that only one of the two mixing 
angles is large. In the case of left-handed squark mixing, for instance, the relevant 
gluino box contribution to the B s -B s mass difference Am s , normalized to the SM 
contribution, is given by [46] 
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where x t = mf/mw, and 



F ( x b L ~ 9 > x hg) + F ( x s 



•Esrg) 2R(xj lr ~, x 



Lgi ^s L g 



b L gi -^SLg, 



C(x) 



X 



12x 3 + 15x 2 — 2x + 6x 3 In x 



(5.11) 
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4(x- l) 3 

For right-handed squark mixing the ratio \Amf R /Am^ M \ would be obtained from 
the above result via the substitution L — > R. 

In Figure 9 the ratio of the SUSY 
contribution to Am s to the SM result is 
shown clS cL function of the mass splitting 
between strange and bottom squarks. The 
same plot with obvious substitutions L — > 
R applies to the case of mixing between 
the right-handed squarks. We observe 
that Am s would greatly exceed the SM 
value in regions of parameter space as- 
sociated with large SUSY contributions 
to the penguin coefficients. To gauge the 
potential impact of a Am s measurement 
near the predicted SM value we impose as 
an example the hypothetical constraint 
that |Amf L /Amf M | < 2. According to 
Figure 9, it then follows that, e.g., | sin 26 l\ < 0.52 for m 
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Figure 9: Ratio \Arrig / Am s \ in units of 
sin 2 2#£ versus the mass splitting m~ s 



m 



in GeV, for m; 



250 GeV. 



I'L 



2mr 



500 GeV, and 



sin 2^1 < 0.28 for m 5 



4mr 



1000 GeV. Hence, if such a constraint would have 



7 L \ \ w.^u iui ,,o SL — -*n %L 

to be imposed in the future, the allowed magnitude of the SUSY contributions to the 
penguin coefficients would be reduced by a significant amount. (Note, however, that 
the coefficients of the chromomagnetic dipole operators are very weakly constrained 
by B s -B s mixing.) 
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Finally, we comment on implications of naturalness for the large right-handed 
up-down squark mass splitting necessary to obtain sizable SUSY contributions to 
the electroweak penguin coefficients. Following [62] we note that in models in which 
SUSY is broken at high energies, e.g., supergravity, there is a naturalness constraint 
on the squark and slepton mass spectrum coming from the hypercharge D-term. In 
these models the Z-boson mass is given by 

m| _ m\ - m 2 2 tan 2 /3 

T ~ tan 2 /? - 1 ' {bA6) 

where m 2 , ml and tan/3 are the usual parameters of the Higgs potential in the 
minimal SUSY extension of the SM [63]. The renormalization-group equations give 



ml = Tr (m% + m\ — 2m? — mf + m~ \ + ml , 

1 oo \ Ql d R ur L l £rJ 2 ' 



— ir [m-Q L +m- iR 



2 \ -1 



where m\ (A = Q L ,u R ,dR, L L ,e R ) are the values of the squark and slepton masses 
at the GUT scale, m 2 , contains the dependence on the other soft SUSY breaking 
masses, and the trace is taken over generation space. It is clear from (5.14) that a 
large mass splitting between the first generation right-handed up and down squarks 
poses a potential naturalness problem. However, the hypercharge D-term can vanish 
in GUT theories in which hypercharge is embedded in the GUT group. For example, 
in SU(5) one has the relations 

2 2 2 2 2 2 2 /crirr\ 

m Q L = m u R = m e R = m l0 , m d R = m L L = m l > I 5 ' 15 ) 

so that it is possible to have large up-down squark-mass splitting without encoun- 
tering difficulties with naturalness. 

To summarize, we have seen that SUSY contributions to the Wilson coefficients 
of the penguin operators in the effective Hamiltonian for b — > sqq transitions, and in 
particular of the isospin-violating electroweak penguin operators, can be substantial 
if the gluino and certain squarks have masses near the weak scale, and other squarks 
have masses near a TeV. Large left-handed or right-handed strange-bottom squark 
mixing is also required. The latter option is naturally realized in certain SUSY 
theories of flavor. 



5.2 Two-Higgs— doublet models 

In extensions of the SM containing charged Higgs bosons [63] there are new pho- 
ton and Z penguin diagrams contributing to the Wilson coefficients of the penguin 
operators. Here we consider a general class of two-Higgs-doublet models (2HDMs) 
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discussed in [64], which contains the conventional type-1 and type-2 2HDMs as spe- 
cial cases. We find that, if terms of order m b /m t are neglected, the new penguin 
contributions only involve the Ht^L coupling. Following [64] we write the corre- 
sponding term in the Lagrangian as 



-Ht R b L 



-it m t b Lj V\ t Ri H + h.c, 



(5.16) 



where V^- is the CKM matrix. In principle, the parameter £ t may contain a CP- 
violating phase. However, the penguin contributions only depend on \(,t\ 2 , and thus 
any weak phase would cancel out. This conclusion holds true even in a wider class of 
multi-Higgs models [65]. It is therefore sufficient to focus on the conventional type-1 
or type-2 two-Higgs-doublet models, for which |£ t | 2 = cot 2 (3. 

It follows from the above discussion that in 2HDMs there are no New Physics 
contributions to the CP- violating parameters p and b entering the parametrization of 
the 5 ± — > nK decay amplitudes in (3.1). Therefore, it is sufficient for our purposes 
to focus on the electroweak penguin coefficients, which induce a new contribution to 
the parameter a. Including both the photon and Z penguin diagrams, one obtains 
at the weak scale [66, 67] 



Co = 



~8T 



where x t n 



8tt 
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In this model there is a simple result for the New Physics contribution to the 
parameter a, normalized to the SM contribution. We find 



- cot 2 p 



EW 



f(xw) 

3 In xt 



1 + 



1 -0.74sin 2 # 



w 



xt — l 



1 + 



9{xm) 
9x t f(x tH ) 



(5.19) 



where the first term in parenthesis is free of hadronic uncertainties, and the second 
term is numerically small. 

In Figure 10, we show the ratio cinp/^ew in units of cot 2 /? as a function of the 
Higgs mass. Even for cot (5 — 1 a significant contribution to the parameter a requires 
a small Higgs mass. This possibility is not excluded by direct searches, however in 
the context of specific models the constraint from the B — > X s 7 branching ratio 
often favors a larger mass and cot 2 /? < 1 [68, 69]. Therefore, it appears unlikely that 
a large new contribution to a can be obtained in 2HDMs. 
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5.3 Models with anomalous gauge- 
boson couplings 

The SU(2)l x U(l) y gauge symmetry of 
the SM fully determines the form of the 
dimension-4 operators that describe the 
vector-boson self-couplings. New Physics 
may induce anomalous couplings of the 
electroweak gauge bosons, which at low 
energies give rise to higher-dimensional 
operators, whose effects are suppressed 
by inverse powers of the New Physics scale 
A. The effects of anomalous gauge-boson 
couplings on rare B decays have been in- 
vestigated in [40, 70], and their impact on 
the determination of 7 from B^ 1 — > nK decays has been discussed in [19] . A general 
parametrization of the anomalous gauge-boson couplings can be found in [71]. In 
low-energy processes, the four new parameters that enter are Ak 7 , Agf , A 7 and g§ . 
The first two represent corrections to couplings already present in the SM, whereas 
the latter two refer to new vertices. Note that these four parameters are real and 
thus can only contribute to the quantity a, but not to p and b. As in the previ- 
ous section, we therefore focus only on the coefficients of the electroweak penguin 
operators. They are [72] 



Figure 10: Ratio onp/^ew m units of 
cot 2 /3 in a general 2HDM, as a function of 
the charged-Higgs mass. 
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with 



f A (x) = -3A<7fln 



h A (x) 



Ak 7 , A 2 
In 




(5.21) 



where the New Physics scale A acts as an ultraviolet cutoff. 

As in the case of the 2HDMs, there is a very simple result for the New Physics 
contribution to the parameter a, normalized to the SM contribution. We find 
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(5.22) 



where as before the first term in parenthesis is free of hadronic uncertainties. For 
instance, taking A = 1 TeV gives 

fl NP _ . a nr\ A Z n AV.„Z 



'EW 



4.20A#f - 0.45# 5 Z + 0.19Ak 7 + 0.03A" 



(5.23) 
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From naive dimensional analysis, one expects that Agf , An' ~ (g w v / A) 2 ~ 10~ 2 
(with v ks 246 GeV the Higgs vacuum expectation value), whereas gf and A 7 are 
expected to be further suppressed [70, 73]. Potentially the most important contri- 
bution in (5.23) is due to Agf, which is bounded by experiment to lie in the range 
—0.113 < Agf < 0.126 [74]. Even when this bound is saturated, Onp/^ew cannot 
exceed 0.5 in magnitude. However, from naive dimensional analysis a^p is naturally 
an order of magnitude smaller. 

6. Conclusions 

We have explored how New Physics could affect purely hadronic FCNC transitions of 
the type b — > sqq focusing, in particular, on isospin violation. Unlike in the Standard 
Model, where isospin-violating effects in these processes are strongly suppressed by 
electroweak gauge couplings or small CKM matrix elements, in many New Physics 
scenarios these effects are not parametrically suppressed relative to isospin-conserving 
FCNC processes. In the language of effective weak Hamiltonians, this implies that 
the Wilson coefficients of QCD and electroweak penguin operators are of a similar 
magnitude. For a large class of New Physics models, we find that the coefficients of 
the electroweak penguin operators are, in fact, due to "trojan" penguins, which are 
neither related to penguin diagrams nor of electroweak origin. 

We have calculated the Wilson coefficients of the penguin operators in the effec- 
tive weak Hamiltonian in several New Physics models, extending the usual operator 
basis where appropriate. We have also included penguin operators mediating the 
decay b — > dsd, which is highly suppressed in the Standard Model. Specifically, we 
have considered: (a) models with tree-level FCNC couplings of the Z boson, extended 
gauge models with an extra Z' boson, SUSY models with broken R-parity; (b) SUSY 
models with R-parity conservation; (c) two-Higgs-doublet models, and models with 
anomalous gauge-boson couplings. In case (a), the resulting electroweak penguin 
coefficients can be much larger than in the Standard Model because they are due to 
tree-level processes. In case (b), these contributions can compete with the Standard 
Model coefficients because they arise from strong-interaction box diagrams, which 
scale relative to the Standard Model like (a s /a)(m^/mg USY ). In models (c), on the 
other hand, isospin-violating New Physics effects are not parametrically enhanced 
and are generally smaller than in the Standard Model. 

We have focused on the rare hadronic decays B ± — > irK, which are particularly 
sensitive to isospin-violating effects. These decays are especially useful for probing 
New Physics contributions, since in the Standard Model the theoretical description 
of such effects is very clean. We have found that the ratio of the CP-averaged 
branching ratios defined in (3.11) and the value of the weak phase 7^ extracted from 

— > ttK decays are the most useful observables for probing isospin-violating New 
Physics contributions. Using a fully general parametrization of the decay amplitudes, 
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we have derived model-independent bounds on in the presence of New Physics, and 
have investigated by how much 7^ could differ from the true value of the CKM angle 
7. We have seen that, depending on the measured value of i?*, it may be possible 
to unambiguously distinguish between isospin-violating and isospin-conserving New 
Physics contributions. Irrespective of the value of R*, we find that large shifts in 7 can 
be caused by even moderate isospin-violating contributions to the decay amplitudes 
of order 10%. In contrast, significant shifts due to isospin-conserving New Physics 
effects would require a new 0(1) CP- violating contribution to the amplitudes. 
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Table 1: Maximal contributions to the relevant phenomenological parameters, as defined 
in the text. Entries for p containing a "— " are unconstrained. Entries marked with a "*" 
are upper bounds obtained using the current values of e 3 / 2 an d R*- For the case of SUSY 
with R-parity the first (second) row corresponds to maximal right-handed (left-handed) 
strange-bottom squark mixing. For the 2HDM we take m#+ > 100 GeV and tan/? > 1. 

For each New Physics model we have explored which regions of parameter space 
can be probed by the B ± — > nK observables, and how big a departure from the Stan- 
dard Model predictions one can expect under realistic circumstances. In Table 1, we 
summarize our estimates of the maximal isospin-violating and isospin-conserving 
contributions to the decay amplitudes, as parameterized by |cznp + ib\ and \p\, re- 
spectively. For comparison, we recall that in the Standard Model a ~ 0.64 and 
b pa p ps 0. We also list the corresponding maximal values of the difference \"f n K 
As noted above, in models with tree-level FCNC couplings New Physics effects can 
be dramatic, whereas in SUSY models with R-parity conservation isospin-violating 
loop effects can be competitive with the Standard Model. In the case of SUSY mod- 
els with R-parity violation, we have derived interesting bounds on combinations of 
the trilinear couplings \' ijk and \'- jk , which are given in (4.24) and (4.25). 

It is worth pointing out that isospin- or, more generally, SU(3) flavor-violating 
New Physics effects in hadronic weak decays could also be important in other pro- 
cesses. For instance, they have been shown to yield a potentially large contribution 
to the quantity e'/e in K — > tttt decays [15]. Moreover, there are other B and B s 
decay channels that could be sensitive to flavor-violating New Physics contributions. 
We look forward to returning to this subject in an extra dimension. 
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